Dynamical regimes and hydrodynamic lift of viscous vesicles under shear 
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The dynamics of two-dimensional viscous vesicles in shear flow, with different fluid viscosities r\ ln 
and ryout inside and outside, respectively, is studied using mesoscale simulation techniques. Besides 
the well-known tank-treading and tumbling motions, an oscillatory swinging motion is observed in 
the simulations for large shear rate. The existence of this swinging motion requires the excitation 
of higher-order undulation modes (beyond elliptical deformations) in two dimensions. Keller-Skalak 
theory is extended to deformable two-dimensional vesicles, such that a dynamical phase diagram can 
be predicted for the reduced shear rate and the viscosity contrast Jjin/^out- The simulation results 
are found to be in good agreement with the theoretical predictions, when thermal fluctuations are 
incorporated in the theory. Moreover, the hydrodynamic lift force, acting on vesicles under shear 
close to a wall, is determined from simulations for various viscosity contrasts. For comparison, the 
lift force is calculated numerically in the absence of thermal fluctuations using the boundary-integral 
method for equal inside and outside viscosities. Both methods show that the dependence of the lift 
force on the distance y crn of the vesicle center of mass from the wall is well described by an effective 
power law y~^ for intermediate distances 0.8R P < y cm < 3-R P with vesicle radius R p . The boundary- 
integral calculation indicates that the lift force decays asymptotically as l/[y C m ln(j/ cm )] far from the 
wall. 



PACS numbers: 87.16.D-, 82.70.-y, 47.15.G- 



I. INTRODUCTION 

Vesicles are fluid droplets enclosed by a fluid lipid 
membrane. Typically, vesicles have sizes of the order 
of 100 nanometers to 10 micrometers, whereas the thick- 
ness of the membrane is only of the order of a nanome- 
ter. Therefore, the membrane can often be regarded as 
a two-dimensional manifold. Vesicle shapes and fluctua- 
tions are then governed by the curvature elasticity. This 
description has been very successful to explain vesicles 
behavior in thermal equilibrium 

The dynamics of fluid vesicles in shear flow has at- 
tracted much attention recentl y [p , d, d, H H 0, H H 
03, HO, El, El EI EE EE El Ej Aspherical vesicles 
under shear can be found in different dynamical phases, 
depending on the viscosities rji n and ?7out of inner and 
outer fluids, respectively, the membrane viscosity ?7 m b, 
the bending rigidity k, the shear rate 7, the membrane 
area, and enclosed volume. As long as shape relaxation 
times of the vesicle are small compared to the time scale 
set by the shear rate 7, the vesicle is always close to its 
equilibrium shape. Under these conditions, vesicles can 
be found either in tank-treading (TT) motion, or - if the 
viscosity contrast A = r]i n /r) out exceeds a critical value - 
in a tumbling (TB) motion. In the tank-treading regime, 
the vesicle shape and orientation are stationary in time, 
but the membrane rotates around the vesicle's center of 
mass in the same direction as the rotational part of the 
shear flow. Here, the orientation is characterized by the 
inclination angle 9 with respect to the flow direction. In 
the tumbling regime, the long axis of the vesicle performs 
a periodic rotation. Keller and Skalak Q developed a 



theory for fluid vesicles with fixed ellipsoidal shape and 
different viscosity contrasts, which is able to explain the 
observed experiments. In recent years, computer simu- 
lations [1 U, H H 0] have shown that the Keller-Skalak 
(KS) theory provides indeed a very good description of 
tank-treading and tumbling. 

However, the vesicle dynamics is far less understood 
when the shear rate is large enough that the vesicle 
cannot relax into its equilibrium shape. Only recently, 
it was shown that a third dynamical regime can ap- 
pear under these conditions, the swinging (SW) regime 
lU E1E1, EI El El, El E3 ~ also called the trem- 
bling [l lj or vacillating-breathing regime [l3||. In the 
swinging state, oscillations of shape and inclination an- 
gle together determine the vesicle dynamics. Swinging 
vesicles were first observed experimentally in Ref. [111 ]. 
With increasing shear rate, a transition from tumbling to 
swinging motion was found. A perturbation theory for 
quasi-spherical vesicles to lowest order in the deviation 
from the spherical shape predicted swinging for a range of 
viscosity contrasts [131 ] : however, since the shear rate ap- 
pears only as basic (inverse) time scale in this approach, 
the experimental results could not be explained. There- 
fore , highe r-order expansions for quasi-spherical vesicles 
0, IH El and a generalized Keller-Skalak (KS) theory 
for ellipsoidal vesicles [l5[ have been developed, which 
are able to predict phase transitions with varying shear 
rate and thereby to explain the experiments of Ref. [ll[ . 

The dynamics in the TT, TB, and SW phases has been 
studied mainly for single vesicles in an unbounded fluid. 
However, in particular due to its physiological impor- 
tance, it is of high interest to study the dynamical be- 
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havior of vesicles under shear in the presence of walls. In 
this case, vesicles are repelled from a wall due to a hydro- 
dynamic lift force Fl . The hydrodynamic lift force plays 
an important role in circulatory systems of vertebrates. 
Since the lift force pushes red blood cells to the center 
of a blood vessel, where the flow velocity is largest, it in- 
creases the efficiency of oxygen transport. On the other 
hand, white blood cells move along the vessel walls in 
order to find defects in the vascular endothelium [IE [2(| • 
This is achieved by special ligands, which are located at 
the outside of white blood cells and bind to receptors 
on the vessel wall to resist the hydrodynamic lift force 




The existence of a hydrodynamic lift force was first 
reported by Poiseuille in 1836 [23|, who observed this 
effect on blood cells. In recent years, the hydrody- 
namic lift force was studied intensively, both theoreti cally 
0, M, EE H3, S] and experimentally 0, HE ED, Hf . 
Abkarian et al. [IE UK observed the unbinding of a heavy 
vesicle, which was pulled by gravity towards a wall, with 
increasing shear rate. For vesicles which are not in direct 
contact with the wall, only studies in three dimensions 
with equal viscosities of inside and outside fluids exist. 
Both, boundary- integral simulations [28| as well as the- 
oretical studies [HI, HH show that the lift force decays 
with a power law 1 /y^ m with increasing distance between 
the vesicle's center of mass and the wall. For vesicles in 
two dimensions, there are only theoretical and numeri- 
cal studies which focus on adhering vesicles bound to the 
wall by a short-ranged attractive potential [IE (2?| ■ 



In this paper, we study the dynamics of a two- 
dimensional (2D) vesicle as a function of viscosity con- 
trast A and shear rate 7, both in the bulk and near a wall. 
The advantage of simulations of a vesicle in two dimen- 
sions is (i) the reduced numerical effort of hydrodynamics 
simulations, which allows for larger system sizes, longer 
accessible time scales, and better statistics, and (ii) the 
simpler form of the equations of the KS theory, where 
no integrals remain in the geometric factors - unlike in 
the 3D version (see App.|A}. This facilitates a detailed 
comparison of the results of theory and simulations. Us- 
ing a mesoscopic hydrodynamics approach, we first show 
that the SW mode also exists in two dimensions, and de- 
termine the dynamical phase diagram. The simulation 
results are compared with the predictions of a general- 
ized KS theory. Second, we study the lift force Fl of 2D 
vesicles, by covering the full range of wall distances y CTa , 
and investigate the effects of viscosity contrast A. More- 
over, we investigate the effect of a wall on the TT-TB 
behavior. For comparison with the results of mesoscopic 
hydrodynamics simulations, we also determine and 
the inclination angle 8 by the boundary-integral method 
for tank-treading vesicles with A = 1. 



II. THEORY AND METHODS 
A. Dimensionless Parameters 

In a 2D vesicle, the perimeter L p and the enclosed area 
A are kept constant (analogously to the constant mem- 
brane surface and the enclosed volume of 3D vesicles) . It 
is useful to combine these two parameters into a dimen- 
sionless quantity, the reduced area 

Here R p = L p /2tt and Ra = \J~ApK are the radii of 
circles with the same L p and A as those of the vesicle, 
respectively. A* is the ratio between the enclosed area 
A and the area of a circle with the same perimeter L p . 
We focus here on a reduced area of A* = 0.7 as a repre- 
sentative for vesicles which deviate significantly from the 
circular shape. 

Shape and orientation of the vesicles are quantified 
by a shape parameter a and inclination angle 9 based 
on the gyration tensor of the vesicle membrane. When 
A max and A m j n are the two eigenvalues of the gyration 
tensor (A max > A min ), and e max and e min the corre- 
sponding eigenvectors, the "asphericity" is described by 
a = (A max - A min )/(A raax + A min ) and the vesicle orien- 
tation by the inclination angle 9 = Z(x, e max ), where x 
is the shear and y the gradient direction. 

The stability of dynamical phases mainly depends on 
two parameters, the viscosity contrast A and the reduced 
shear rate 

;= TWg (2) 

K 

The time rj ou tR p /n is the characteristic relaxation time 
in thermal equilibrium, where k is the bending rigidity. 
Thus 7* expresses the interplay between the perturbation 
by the external field 7 and the ability of the vesicle to 
restore its equilibrium shape. 



B. Generalized Keller-Skalak Theory in Two 
Dimensions 

Keller-Skalak (KS) theory is based on the assump- 
tion that vesicles have a fixed ellipsoidal shape. There- 
fore, it cannot describe the swinging state with oscillating 
vesicle shapes. Therefore, KS theory has been general- 
ized to include shape deformation in three dimensions 
[lH . This theory is applicable to ellipsoidal vesicles over 
a wide range of reduced volumes, while higher-order per- 
turbation theory [IE E3, EH ^ s limited to quasi-spherical 
vesicles. Here, we employ the two-dimensional version of 
the generalized KS theory. The differential equations for 
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the asphericity a and inclination angle 9 are given by 
1 da 
~^~dt 



An area potential 



b Q R p dF . 

1 Disin(20), (3) 



A*7 



dc 



~ = ^[-1 + 5(0)008(2(9)], 



with prcfactors 



and b\ = 



(4) 



(5) 



u 4tt(A + 1) 1 2(A + 1) 

There are no adjustable parameters. An explicit expres- 
sion for B(a) and its derivation are described in App. [Al 

The time evolution of 9 is described by Eq. (j4|), which 
has the same form as in two-dimensional KS theory. 
However, B{a) is now not constant but depends on the 
time-dependent vesicle shape a(t). The time evolution 
of a (see Eq. (J3J) is derived based on the perturbation 
theory of quasi-circular vesicles [33|. Here, F is the free 
energy of the vesicle shape at constant A* . F attains its 
minimum for an elliptical vesicle shape in equilibrium. 
Thus, the first term on the right hand side of Eq. ([3]) 
causes a relaxation of a towards its equilibrium value. 
The second term represents the change of a due to the 
external flow field. Eqs. (|3|) and (j4|) arc solved numeri- 
cally using a fourth-order Runge-Kutta method. 

Thermal fluctuations can be incorporated in this ap- 
proach by adding Gaussian white noises g a {t) and ge(t) 
to Eqs. ^ and respectively The noise terms obey 
the fluctuation-dissipation theorem, such that (gi(t)) = 
and {g l {t)g 1 {t')) = (2k B T/Q)6 tJ 5{t - t') with i,j G 
{a, 9}, where k B T is the thermal energy. As a reasonable 
approximation, we employ the rotational friction coeffi- 
cients of a circle, 

47T 

Co = ^-?7outi?A 2 (A + 1) and Q = 47rry out i? p 2 . (6) 



C. Mesoscale Hydrodynamics Simulation Method 

1. Membrane Model 

The membrane is modeled by a closed chain of n 
monomers of mass M. For a monomer with index i (with 
1 < i < n), we introduce the notation 

i- = (i — 1) mod n and i+ = (i + 1) mod n (7) 

for the indices of its two neighboring monomers. Thereby, 
the ring topology is taken into account correctly. The 
monomers are connected by a harmonic spring potential 



h 



i=l 



(8) 



where := r; + — are the bond vectors, and I is the re- 
laxed bond length. The curvature elasticity of the mem- 
brane is described by the bending potential 



(9) 



U A = ^(A-A )\ 



(10) 



is introduced to control the deviations of the area A from 
its target value A . Here, the enclosed area A in Eq. (fTO|) 
is obtained from the monomer positions by 
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n 



(11) 



i=l 



2. Multi-Particle Collision Dynamics 

For the solvent hydrodynamics, we employ multi- 
particle collision dynamics (MPC), a particle-based 
mesoscopic simulation technique HH, [H, |3(|. The dy- 
namics of an MPC fluid evolves in two alternating steps. 
In the "streaming step" , particles move ballistically for a 
time At, the collision time, according to their current ve- 
locities. For the "collision step" , solvent particles are first 
sorted into the cells of linear size a of a regular square lat- 
tice; all particles in a cell then exchange momenta such 
that the total translational momentum is conserved in 
each collision cell. 

Several modifications of the orig inal MPC algorithm 
have been introduced recently [37| . which differ in the 
way the collision step is executed. We employ the MPC- 
AT+a version of multi-particle collision dynamics, which 
uses an Anderson thermostat (AT) and locally conserves 
angular momentum (+a) in addition to translational mo- 
mentum. In MPC- AT, new particle velocities relative to 
the center-of-mass velocity are chosen from a Maxwell- 
Boltzmann distribution with temperature T. This ther- 
mostat avoids any heating due to energy dissipation in 
sheared system. For details of the MPC-AT+a algo- 
rithm, see Refs. [13, H^. We use this algorithm, since 
local angular-momentum conservation is crucial in binary 
fluid systems with different viscosities [3!| . 

Simulations are performed with a rectangular simula- 
tion box with linear sizes L x and L y , periodic boundary 
conditions in the x direction, and no-slip wall bound- 
ary conditions in the y direction. Linear shear flow with 
shear rate 7 is realized by moving the upper wall with a 
velocity jL y x, whereas the lower wall is held at rest. 

Many properties of the MPC-AT+a solvent can be ad- 
justed by the simulation parameters collision time At, 
the particle number density n s , and the particle mass m. 
The solvent viscosity r\ = rjkin + Vcott is a sum of a kinetic 
rykin and a collisional contribution 7j co ii, which have been 
calculated analytically (38j . 



Vkiz 



Vcol 



n s k^TAt 



m(n s 



7/5) 



24At 



(12) 



(13) 



The viscosity r] out of the fluid outside of the vesicle is ad- 
justed by varying the collision time At in the range from 
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At = 0.003a vW fc BT to 
these collision times the 
than the cell size a, the 
nated by 77 co n (see Eqs. 
sional viscosity ?7 co n oc m 
varied by using different 
and outer fluid particles, 

A = 



At = 0.01a^m/fc B T. Since for 
mean free path is much smaller 
total shear viscosity r\ is domi- 
(fT2"j) and (fT5jO, Since the colli- 
the viscosity contrast A can be 
masses mj n and m of the inner 
respectively, which implies 



Vout 



m 
m 



in 



(14) 



In our simulations, the viscosity contrast is varied from 
A = 1 to A = 10 (with m < m ln < 10m), while all the 
other MPC parameters are the same for the fluid on both 
sides of the membrane. 



3. Membrane Interactions 

In order to describe an impermeable membrane in flow, 
it has to be ensured that MPC particles stay on the cor- 
rect side of the membrane (i.e. inside or outside of the 
vesicle). For numerical efficiency, it is advantageous to 
relax this condition for short length and time scales, as 
it was done in previous 3D vesicle simulations @. The 
streaming and collision steps for the fluid particles are 
carried as in the absence of the membrane. This implies 
that after each streaming step, some MPC particles have 
crossed the membrane. For the (few) particles which are 
now located on the wrong side of the membrane, with a 
direction of their velocity which would bring them away 
even further away from the membrane, the velocities have 
to be modified such that they move towards the mem- 
brane instead, in order to cross back to their correct side. 
We denote this velocity update a "membrane collision". 
It has to be constructed such that the translational and 
angular momentum as well as the kinetic energy of the 
fluid particles and membrane monomers are conserved 
locally. Our procedure for membrane collisions is a gen- 
eralization of the standard bounce-back rule for no-slip 
boundary conditions. A detailed description of this pro- 
cedure is provided in App. [B] 

In order to prevent the membrane from crossing the 
walls, a purely repulsive Lennard- Jones potential 




e, 0<y< y^cr 
otherwise 



is employed, which depends only on the distance y of a 
monomer from a wall. 

For the determination of hydrodynamic lift forces, we 
employ a gravitational body force fc = —ygAg, which 
acts on the internal fluid of the vesicle. Here, g denotes 
the strength of the gravitational field, and Ag is the mass- 
density difference between the inner and outer fluids. The 
gravitational body force fo acting on the inner fluid can 
be expressed as a potential Uq, which only depends on 



the monomer positions, 
U G 



6A 4^ 



Vi+Vi + ) 



(15) 



Here, and are the y components of the monomer 
positions and r^ + , respectively, and Fg — \J A fodA\ is 
the total gravitational force acting on the vesicle. Fq has 
a constant value and is used as a simulation parameter. 

As long as not specified otherwise, the parameters used 
in our vesicle simulations are n = 50, I = a = \/2<7, 
n s = 10a~ 2 , M = 10m, e = 10fc B T, and k/1 = 50fc B T. 
For the reduced area, we require that it deviates less 
than 1% from its target value of A* = 0.7. Since A* is a 
function of the perimeter L p and the enclosed area A (see 
Eq. ([1])), the parameters k sp and &a for the potentials 
U sp and Ua, respectively, have to be sufficiently large. 
We chose k sp = 10 4 k B T/a 2 and k A = 80k B T/a 4 . With 
these parameters, the effective vesicle radius is obtained 
to be R p = 7.8 1. The size of the simulation box is L x = 
L y = 80a. Gravitational forces Fq are only applied in 
simulations for the hydrodynamic lift force, where values 
in the range k B T/a < Fq < b0k B T/a are investigated. 

In simulations, different reduced shear rates 7* can be 
achieved, according to Eq. @, by varying 7, 7y out , R pi 
or n. Since equilibrium properties like the undulation 
spectrum depend on R p and K, we vary 7* by adjusting 
7 and 7y ou t. In order to avoid inertial effects, we restrict 
the shear rates to obtain low Reynolds numbers Re = 
7pi?p/77 ou t, where p is the density of the outer fluid. The 
maximum Reynolds number is Re= 0.17. 



D. Boundary-Integral Method 

For comparison with our MPC simulation results of the 
lift force, we also perform numerical boundary-integral 
calculations. The hydrodynamic lift force in 2D has been 
studied previously with the boundary-inte gral approach 
for vesicles in direct contact with the wall [2fl (27[ ■ This 
method has the advantage that it can be used to cal- 
culate lift forces on vesicles even for very large distances 
y cm from the wall and for reduced areas A* close to unity, 
which are not easily accessible by MPC simulations. On 
the other hand, our boundary-integral calculation is re- 
stricted to elliptical shapes and ignores thermal fluctu- 
ations, which give rise, e.g., to undulation-induced re- 
pulsion near a wall. We focus on tank-treading elliptical 
vesicles without viscosity contrast, i.e. A = 1. In the 
steady tank-treading state the lift force can calculated 
from a single, time-independent vesicle shape. Whereas 
in MPC simulations, the wall distance y cm is calculated 
for a given strength of the gravitational force, we fol- 
low the opposite procedure with the boundary-integral 
approach, by calculating the lift force for a given wall 
distance y cm . 

For ellipse half axes a\ and 02, wall distance y cm _, and 
inclination angle 8, the location r of the vesicle membrane 
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FIG. 1: For the boundary-integral calculation, the vesicle is 
discretized into segments i of length U , which have the orien- 
tation xii and a center-of-mass position n. 



is uniquely defined (with x c 
zation by the angle ip, it is 



0). With a parameteri- 



r(p) = y cm y + 



cost 
smt 



— smt 
cos 9 



r'(<p) (16) 



The difficulty is that the force distribution f(r') along 
the membrane is a priori unknown. Instead, we know the 
velocities v(r) at each site of the membrane. Eq. I]20p 
is thereby a Fredholm integral equation of the first type. 
This integral equation is solved numerically. For this pur- 
pose, we discretize the membrane in N straight segments, 
which have to be small enough such that the difference 
in velocities between two neighboring segments is small 
and the force distribution can be assumed to be constant 
along the segment. A segment with index i has a velocity 
Vi, center-of-mass position r i; length and orientation 
Ui (see Fig. [TJ. The discretized form of Eq. ([20)1 is 



2^ / G{r i ,r j +u j s)f(T j + Vi j s)ds. (21) 



J=l -«J 



Since the force distribution fj = f(r^ + Ujs) is assumed 
to be constant over the whole segment j, it can be moved 
outside of the integral, 



where r'(tp) is the membrane position in the principal- 
axis system of the ellipse, 



a\ cos tp 
ai sin ip 



(17) 



In the steady tank-treading state, the center-of-mass 
velocity v cm = xi> cm of the vesicle has only a non- 
vanishing component in shear direction. If v cm and the 
tank-treading angular velocity lj arc known, the velocity 
of the tank-treading membrane is given by 



v(tp) = a 
with the tangent vector 



R p Loi(p). 



t(y>) 



ai 

_ a-2 



— cos 6 

0,2 
0.2 



\r'{p)[ 



(18) 



(19) 



In a tank-treading membrane in shear flow, forces arising 
from pressure and viscous stress have to be balanced in 
order to maintain a steady motion. The force distribution 
f(r') along the membrane dA is related to the velocity 
field at position r by 



v(r) - 7yx 



G(r, r'(s))f{r'(s))ds 



(20) 



■c). 1 



Here, ds is a line element of the membrane dA, and the 
second-order tensor G(r, r'(s)) is the Greens function of 
the Stokes equation which satisfies the boundary condi- 
tions. For vesicles in an unbounded fluid, G(r, r'(s)) is 
the Oseen tensor. In our case of a vesicle near a wall, the 
half-space Oseen tensor - also known as Blake tensor - 
is convenient, as it realizes no-slip boundary conditions 
at the wall. The full expression for the two-dimensional 
Blake tensor can be found in Ref. 14011. 



Vi-72/iX 



N 

E 



W2 



G(r i; Tj + \ijs)ds 



N 



f 



(22) 

The calculation of H^- can be performed analytically, 
both for the free-space and the half-space Oseen tensor. 
Thus, the integral equation ([20)1 is reduced to a set of 
linear algebraic equations which can be easily solved nu- 
merically. 

The segment velocities Vj depend linearly on lj and 
v cm (see Eq. (JTSJl). Therefore, we can extend the linear 
system of equations (|20j) by two additional conditions, 
which determine u and v cm self-consistently in the steady 
state. For the first condition, we require that the sum 
of tangential forces along the membrane vanishes. The 
second condition is that the vesicle does not experience 
a net force in shear direction. The total system of linear 
equations finally reads 



Hyfj - V cm X - tjijp 

3 = 1 



N 



y^x ■ 

i=l 
N 



(23) 
(24) 
(25) 



This set of equations is solved numerically with up to 
N = 600 segments. Once, u>, w cm and the force distri- 
bution are known, quantities like the lift force Fl and 
the torque M on the vesicle, as well as the velocity v(r) 
and pressure fields p(r) in the surrounding fluid can be 
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calculated as 



iV 



3 = 1 
N 



p(r) 



j'=i 

TV 

j'=i 



JV 

E 

i=i 



/ G(r, r_, + Ujs)di 

J-lj/2 

em* 

/ g(r, Tj + u J s)ds 
J-U/2 



(26) 
(27) 

f„(28) 
(29) 



Here, g(r, r') is the half-space pressure vector (sec 
Ref. [13]). The lift force F L and the torque M (sec 
Eqs. Q26p and l|2?p). are thereby functions of the four pa- 
rameters 9, yem, o,\, and a2, which define the membrane 
location uniquely. Using a numerical root finder (Brent's 
method), the stable inclination angle 9, for which the 
torque vanishes, is determined while keeping the other 
parameters y cm , ai, and 0,2 fixed. 



III. DYNAMICAL REGIMES OF VISCOUS 
VESICLES IN UNBOUNDED SHEAR FLOW 

A. Phase Diagram 

We consider first the dynamics of vesicles in shear flow, 
far from walls and in the absence of a gravitational field. 
The 2D generalized KS theory predicts a phase diagram, 
see Fig. [3J which shows the qualitatively the same fea- 
tures as a function of 7* and A as the 3D version [l5[ . At 
small and large A, a vesicle exhibits tank-treading (TT) 
and tumbling (TB) motion, respectively. At large 7* and 
intermediate A, the swinging (SW) phase appears. As in 
3D generalized KS theory, TT with negative inclination 
angles 9 < appears close to the TT-SW transition line. 
The coexistence of two TT states (one with 9 < and 
the other with 9 > 0) or of a TT and a SW states are 
also seen. In 2D, TT with 9 < is stable, unlike in 3D 
[IH, where the vesicle can escape by turning its longest 
axis into the vorticity direction. 

MPC simulation results of the three dynamical regimes 
(TT, TB, SW) are illustrated by a sequence of snap- 
shots in Fig. [3l Fig. [3ja) shows a tank-treading vesi- 
cle, which has a constant shape and orientation (except 
for its thermal membrane undulations). However, as a 
marker on the membrane indicates, the membrane ro- 
tates around the center of mass. A tumbling vesicle is 
shown in Fig. G3Jb) , where the shape remains almost un- 
changed, but the orientation steadily rotates. The mo- 
tion of a marker on the membrane shows that the mem- 
brane is not completely fixed with respect to the vesicle 
shape. Fig. [3Jc) illustrates the swinging state (see also 
movie [HI). As long as the vesicle orientation has a pos- 
itive inclination angle, the elongational part of the shear 
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FIG. 2: (Color online) Phase diagram at A* — 0.7. Dashed 
(red) and solid (blue) lines represent the results of the gener- 
alized Keller-Skalak theory with 61 = 3/[2(A+l)] (see Eq. ©) 
and 61 = 1/(A + 1), respectively. The black circles (•) indicate 
the location of the simulation which are shown in Fig. [5] 



flow causes an elongation of the vesicle shape (increasing 
a). For a large shape parameter a, the vesicle is tem- 
porarily in the tumbling regime, until a negative inclina- 
tion angle 9 is reached (for < tj < 10 in Fig.^c)). For 
negative 9, the elongational component of the flow acts 
to reduce a. Due to the constraint of fixed perimeter L p 
and fixed enclosed area A, the vesicle assumes a potato- 
like shape, such that a decreases (for 10 < tj < 15 in 
Fig. GUJc)). The vesicle is then stretched again by the 
elongational flow leading to a positive inclination angle 9 
and increasing shape parameters a (for tj > 20). 

It is important to note that elliptical deformations are 
not sufficient in 2D, because the constraints on perimeter 
and area complete determine the elliptical shape [33[ - 
in contrast to 3D, where the deformation in the vortic- 
ity direction provides sufficient degrees of freedom [l5[. 
Thus, higher-order undulation modes beyond elliptical 
deformation arc required, which can be seen clearly in 
Fig. EDJc). This is reminiscent of the behavior of 3D vesi- 
cles in an elongational flow after flow reversal [12, l43| . 
where also higher-order undulation modes play an im- 
portant role. 



B. TT-TB Transition 

The generalized KS theory predicts that for small shear 
rates, with 7* < 6, the TT-TB transition at A ~ 3.25 
hardly depends on 7* (see Fig. [2]). In this regime, shape 
deformations are very small, and the behavior can be well 
described by t he original KS theory. We choose a shear 
rate 7 = 0. 01 y/k^T/ ma 2 in our simulation, correspond- 
ing to a small reduced shear rate 7* = 3.6. 

In Fig. [U the dependence of the average inclination an- 
gle 9 on the viscosity contrast A is shown. Our MPC sim- 
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(b) Swinging: 



(c) Tumbling: 
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FIG. 3: (Color online) Sequences of vesicle snapshots for each 
of the dynamical regimes, shown (a: TT, b: SW, c: TB). 
A (red) bullet marks one fixed membrane element to indi- 
cate the membrane motion. All systems share the param- 
eters k/1 = 50/cbT and A* = 0.7. Further parameters are 
(a) ?7out = 36\/kB.Tm/a, A 
(b) 



?7out = SQ^Jk^Tm/a, 
responding to 7* = 3.6; 



1, 7 

»7out 



cor- 
4, 



120V fc B Tm/a, A 
7 = 0. 0333 \Jk,B,T /m/a corresponding to 7* = 38; and ( 
?7out = 3Q\/k-BTm/a, A = 10, 7 = Q.Qly/k&T/m/a corr 
sponding to 7* = 3.6. 
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FIG. 4: (Color online) Inclination angle 6 as a function of vis- 
cosity contrast A for simulations with 7* = 3.6. For compari- 
son the results of the boundary-integral calculation of Beau- 
court et al. [J (without thermal fluctuations) as well as the 
curve of KS theory @] (see Eq. (|30j)) are shown. 



ulations well reproduce the results of previous boundary- 
integral calculations by Beaucourt et al. 0]. Deviations 
close to the TT-TB transition at A* ~ 4 arise from ther- 
mal membrane undulations, which arc present in our sim- 
ulations, whereas the results of Ref. [H have been ob- 
tained in the zero-temperature limit. 

Moreover, thermal fluctuations lead to a continuous 
crossover rather than a sharp TT-TB transition. Thus, 
there arc a few tumbling events already for viscosity con- 



trast A = 3, and also simulations with A > A* ~ 4 exhibit 
some time intervals of tank-treading motion. Our simu- 
lations also show that the existence of a tumbling regime 
depends sensitively on the Reynolds number Re. For Re 
> 1, 6 decreases more gradually with increasing A, and 
no tumbling motion was observed at viscosity contrasts 
as large as A = 10. Thus, we conclude that inertial effects 
enhance the TT-membranc rotation. 

Fig. H] also shows that KS theory Q provides a good 
description of the A dependence of 9 and the TT-TB 
transition. This transition is explained by the KS theory 
as follows. The stationary inclination angle 9 in the tank- 
treading regime is determined by Eq. Q with fixed a as 



1 ( 1 

-arccos I - — 



(30) 



For small A, the inclination angle 9 decreases mono- 
tonically up to a critical viscosity contrast A*, where 
9 = 0. For larger viscosity contrasts A > A*, the tum- 
bling regime, there is no real solution of Eq. (|30[) . i.e. no 
stationary inclination angle exists, and the vesicle per- 
manently rotates. 



C. TB-SW Transition 

To investigate the TB-SW transition, we consider a 
fixed viscosity contrast of A = 4, and perform simulations 
for four different reduced shear rates 7* = 5, 17, 22, and 
38. The locations of these four shear rates in the dynam- 
ical phase diagram are indicated in Fig. [21 The resulting 
trajectories in the 9-a plane arc shown on the left-hand 
side of Fig. [5l In this representation, closed cycles indi- 
cate swinging events, whereas trajectories spanning the 
full [— 7r/2, +7t/2] range of 9 are tumbling events. Obvi- 
ously, thermal noise has a large impact on the vesicle dy- 
namics. In particular, at small inclination angles 9 ~ 0, 
small thermal fluctuations are decisive for the vesicle to 
perform a tumbling or swinging cycle. 

The simulation data of Fig. [5] suggest that the TB-SW 
transition point is located between 7* = 17 and 22. This 
is about a factor 2 larger than the prediction 7* = 9 of 
generalized KS theory. In the generalized KS theory, a 
possible source of error can be found in the estimate of 

60 and 61, which have both been calculated in the cir- 
cular limit, see Eq. (JSJ. Therefore, we also calculate the 
phase diagram with b\ reduced by a factor 2/3, i.e. with 

61 = 1/(A + 1), which gives a better agreement with our 
simulations for non-circular vesicles with A* = 0.7 (see 
Fig. [2]). In this case, the effect of thermal noise on tra- 
jectories is found to be very similar as in the simulations 
for all four reduced shear rates (see Fig. [5]). Therefore, 
the deviations between the generalized KS theory and 
simulations can be alleviated by a small modification of 
prefactors. Further theoretical developments are needed 
to determine the prefactors 60 and b\ analytically for non- 
circular shapes. We conclude that generalized KS theory 



8 



Simulations: Theory: 
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FIG. 5: (Color online) Trajectories in the 0-a plane for A = 
4, both with (thin red lines) and without (thick blue lines) 
thermal noise. The trajectories in (a), (c), (d), and (g) are 
obtained from simulations, whereas the curves in (b), (d), (f), 
and (h) are calculated from the generalized KS theory with 
noise and fei = 1/(A + 1). The reduced shear rates are 7* = 5 
for (a) and (b), 7* = 17 for (c) and (d), 7* = 22 for (e) and (f), 
and 7* = 38 for (g) and (h). In all plots, the corresponding 
theoretical trajectory according the generalized Keller-Skalak 
theory without thermal noise is shown as a thick blue (dark) 
line. 



provides a good description of vesicle dynamics in shear 
flow in both two and three spatial dimensions. 

Elastic capsules [13, 5H and red blood cells [46| can 
also exhibit a swinging motion. However, the angle 9(t) is 
always positive during these oscillations — unlike SW of 
fluid vesicles. The physical mechanism is an energy bar- 
rier for the TT rotation caused by the membrane shear 
elasticity and the anisotropic shape of the spectrin net- 
work O, S3, E3- Although, a vesicle in 2D (a closed 
string) does not have membrane shear elasticity, an en- 
ergy barrier for the TT rotation can be introduced by in- 




FIG. 6: (Color) (a) Shape and (b) pressure field of a tank- 
treading vesicle under shear flow close to a wall, in steady 
state with viscosity contrast A = 3. The color code is ex- 
pressed in units of n s kBT. The hydrodynamic lift force is 
balanced by an external gravitational force Fg = lAk^T/a, 
where the average distance from the wall is j/ cm = 7.96a. See 
also movie [50( ]. 

homogeneities in the spontaneous curvature (49j . In the 
future, it will be interesting to investigate the coupling of 
different swinging mechanisms in composite membranes. 



IV. LIFT FORCE 

We now consider a vesicle under the combined effect of 
a shear flow and a gravitational force Fq , see Fig. EJa) . 
The vesicle moves towards or away from the wall un- 
til gravitational Fq and lift forces Fi,(y cm ) balance each 
other (see also movie |5(j). In this steady state, the lift 
force Fi J (y cm ) equals the gravitational force in magni- 
tude. 

Fig. [BJb) shows the pressure field in the outer fluid for 
the steady-state configuration of a tank-treading vesicle. 
The hydrodynamic lift force is the integral of the pressure 
forces over the membrane contour. The higher pressure 
in the gap between the vesicle and the wall is responsi- 
ble for the lift force. Fig. Efb) also nicely demonstrates 
that there is a lower pressure at the two caps of the vesi- 
cle, which is the origin of vesicle elongation. The hydro- 
dynamic lift force is a pressure force which is of purely 
viscous nature - in contrast to e.g. aerodynamic forces 
acting on the wings of an airplane, which arc caused by 
inertial forces. 

The dependence of the hydrodynamic lift force on the 
wall distance is shown in Fig. [T^a) — calculated as (y C m) 
for fixed gravitational force in the simulations, and as 
lift force at fixed y CTa in the Oseen calculations, as ex- 
plained in detail in Sec. Ill C 31 and Sec. Ill Dl above. Lift 
forces of vesicles with A < 4 can be well described by a 
power-law F^ oc for F^ < 2.5fceT '/ 'R p , correspond- 
ing to y cm > R p . For these distances, the vesicle is not 
in direct contact with the wall. At applied gravitational 
forces larger than 2.5/cb? 1 / R p , the vesicle touches the wall 
(l/cm ^ Rp)- However, the distance y cm between the cen- 
ter of mass and the wall can be reduced even further 
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FIG. 7: (Color online) (a) Lift forces and (b) average incli- 
nation angle as a function of the average wall distance ?/ cm of 
the vesicle, from MPC simulations with 7* = 3.6 and from 
boundary-integral (Oseen) calculations, as indicated. The leg- 
end in (a) applies to both plots. Simulation points marked by 
big (red) circles refer to tumbling vesicles. For comparison, a 
line with the power-law dependence is plotted in (a). The 
right-most data points in (b) correspond to Fg = 0, so that 
they cannot be shown in the double-logarithmic presentation 
in (a). 



by vesicle deformation. The l/y^ m dependence does not 
apply in this regime. Finally, the constraints of fixed 
enclosed area A and fixed perimeter L p keep the wall 
distance larger than y cm > 0.628i? D . 

Fig. E|Ja) shows that the lift forces decrease with in- 
creasing viscosity ratio A for a fixed wall distance y cm . 
This behavior is analyzed in more detail in Fig. where 
the amplitude Fi,y^ m of the lift force is plotted as a 
function of the viscosity contrast A. Although solid col- 
loidal particles of elliptical shape experience no net lift 
force [5l[ , tumbling vesicles with finite A obtain lift force 
due to an asymmetry of its shape deformations and a 
small tank-treading component (compare Fig. |3fc)). 

For vesicles in three dimensions, both boundary- 
integral simulations [28| as well as theoretical studies 
[24l |25|] show a l/y 2 m dependence of the lift force for 
vesicles far from the wall. The theory of Olla 0, [H| as- 




FIG. 8: (Color online) Amplitude of the lift force, 
FL,y% m /(kBTR p ), as a function of the viscosity contrast A. 
The amplitudes are fits to the curves in Fig. for which 
vesicles are not in direct contact with the wall. 



sumes an ellipsoidal shape for the vesicles with half axes 
a-i, a-2, a-3 <C 2/cm- It is not possible in this case to de- 
rive expressions for two dimensions by taking the limit 
as — ► 00 — as done in App. [X] for KS theory of vesi- 
cles in unbounded flows — because this limit is incon- 
sistent with the assumption ai,a2,&3 *C y C m- Therefore, 
instead of an analytical theory, we perform boundary- 
integral calculations of 2D elliptical vesicles with A = 1 
in the presence of a wall, as described in Sec. Ill Dl For 
the results in Fig. Ela), the effect of the opposite wall 
at L y = lORp is also taken into account by plotting 
Fl (j/cm) - Fl (L y - y cm ), where F L (y cm ) and F L (L y ~ y cm ) 
are obtained from two independent boundary-integral 
calculations. Of course, a more precise calculation would 
require the use the two-wall Oseen tensor |4£j instead of 
the half-space Greens function. However, as long as the 
distance L y between the two walls is sufficiently large, 
Fh(y C m) — Fl,(L v — y C m) is a good approximation of the 
two-wall lift force. Fig. [TJa) as well as Fig. [8] show that 
the boundary-integral calculation indeed agrees nicely 
with the corresponding MPC simulation for A = 1. The 
amplitudes Fi,y cm 2 only differ by about 25%. Reasons for 
this deviation are that the boundary-integral calculation 
is done for elliptical shapes, whereas vesicles in simula- 
tions are closer to the equilibrium shape (compare Fig. [T] 
and Fig. [6|). Moreover thermal fluctuations in the MPC 
simulations may cause differences. 

Fig. [7ja) shows that tumbling is suppressed when the 
gravitational force exceeds a threshold value, depending 
on the viscosity contrast A. In order to perform a tum- 
bling motion, the center-of-mass distance y cm has to be 
on the order of or larger than the long vesicle axis a\. 
However, for larger gravitational forces, the center-of- 
mass distance y cra becomes smaller than a\, such that 
even vesicles with high viscosity contrasts do not tumble. 
Even if y cra is slightly larger than ai, the vesicle cap has 
to come so close to the wall that the resulting pressure 
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FIG. 9: (Color online) Dynamic phase diagram of tank- 
treading and tumbling states as a function of the gravita- 
tional force Fg and the viscosity contrast A. The reduced 
shear rate is 7* = 3.6. Circles (•) indicate tumbling motion 
(TB), squares (■) tank-treading motion (TT). The line for 
the phase boundary is a guide to the eye. 
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FIG. 10: (Color online) Hydrodynamic lift force Fl on vesi- 
cles with varying A* vs. wall distance y cul obtained from 
boundary-integral calculations. 



forces prevent the inclination angle to reach tt/2. This 
leads to the dynamical phase diagram shown in Fig. [9] 
When the gravitational force is large, tumbling is sup- 
pressed and the vesicle displays a tank-treading motion 
at the wall. With increasing A, the gravitational force 
necessary to prevent tumbling increases. 

The dependence of the inclination angle 9 on the wall 
distance y cm is shown in Fig. [7]Jb). Without a gravita- 
tional force, the lift force caused by the upper wall at 
y = 10R P compensates the lift force of the lower wall 
when the vesicle is in the center, at y cm = 5R p . Since 
the lift forces are very small nearby, strong fluctuations 
are observed in the wall distances for small Fq . 

As long as a vesicle is tank-treading, its inclination an- 
gle 9 decreases when it approaches the wall. Even if the 
vesicle does not touch the wall, the pressure at the lowest 
part of the membrane is highest (see Fig.[6jb)) such that 
it causes a torque which lowers 9. For very small wall 
distances, the vesicle comes into direct contact with the 
wall, where the repulsive wall potential causes an addi- 
tional torque, which decreases 9 even further, until the 
vesicle is finally completely parallel to the wall. 

Vesicles with A > 3 start to tumble at sufficiently large 
wall distances. Since vesicles with viscosity ratios A = 3 
and A = 4 are still tank-treading most of the time and 
only occasionally perform a tumbling motion, their in- 
clination angles arc non-zero, whereas for A — 10 and 
Fq < 1, the average inclination angle 9 essentially van- 
ishes (see Fig.JTJb)). 

We employ the boundary-integral approach to calcu- 
late the dependence of the hydrodynamic lift force of 
vesicles with A = 1 on the reduced area A* in the ab- 
sence of thermal fluctuations. Also, since this method is 
not restricted by the system size (as simulations) , the lift 
forces can be calculated even for very large wall distances. 
Fig. [TU] shows the hydrodynamic lift force as a function 
of the wall distance y cm for different reduced areas A* in 
the range 0.65 < A* < 0.99. This plot shows that the 
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FIG. 11: (Color online) Amplitude of the lift force Fl as a 
function of 1 — A* . The data points are fits to the curves in 
Fig. [10] for y cnl > 60i? p . These data points (Oseen) are fitted 
to a Vl — A* dependence for small values of 1 — A* . 



y~ni dependence of the hydrodynamic lift force Fl only 
holds for distances y cm < 5R P from the wall. For larger 
distances y cm , a crossover to a 1/ [(y cm /R p ) In (y cm /R p )] 
dependence is obtained. This power law with a loga- 
rithmic correction fits perfectly the numerical data of 
the boundary-integral calculation for y cm > 2R p for all 
considered reduced areas. Fig. [TT1 shows the amplitudes 
K = FLy cm ln(y cm /R p ) of the lift forces Fl vs. 1 — ^4* in 
the far- field limit. These amplitudes K are determined 
by fitting the expression F L = K/ [(y cm /R p ) In (y cm /R P )] 
to all data with y cm > 60F p . For 1 — A* = 0, i.e. for 
a circular vesicle, the lift force vanishes, which directly 
follows from the time-inversion symmetry of the Stokes 
equation. For small deviations from the circular shape, 
the lift force rapidly increases with 1 — A*, and follows 
a power-law dependence K ~ s/l — A* ~ yl A p /7r, with 
the excess length A p = L p /Ra - 2tt, for 1 — A* < 0.1 
(see Fig. [TTj) . We observe a monotonic increase of the 
lift force with increasing 1 — A* over the whole range of 
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reduced areas A*. 

Since boundary-integral calculations do not take into 
account thermal fluctuations, lift forces can be deter- 
mined for large y cm (only limited by numerical accuracy). 
However in simulations as well as in real systems, lift 
forces at large y cnl have a vanishing effect compared to 
thermal noise. Moreover, with very large wall distance 
y cm ~ yj 77 ou t / pj, inertial effects are not negligible, so 
that the Stokes approximation becomes less reliable [2 51 ] . 



V. SUMMARY AND CONCLUSION 

We have studied the dynamics of vesicles in shear flow 
in a two-dimensional model system. This system shows 
a variety of interesting dynamical phenomena. 

First, we have investigated the effect of the viscosity 
contrast A, i.e. the ratio between the inner and outer vis- 
cosities of a vesicle, on the dynamics in unbounded flows. 
With increasing A, the sequence from "tank-treading" 
over "swinging" to "tumbling" motion is generically ob- 
served — except for small shear rates j, where the in- 
termediate swinging phase is absent. Thus, the swinging 
phase appears in the phase diagram of 2D vesicles un- 
der shear in the same way as it was found previously 
for 3D ellipsoidal vesicles. However, the mechanism of 
swinging is different in two and three dimensions. While 
in 3D, ellipsoidal deformations are sufficient to obtain 
swinging, in 2D higher-order undulation modes are re- 
quired. Thermal fluctuations play an important role; 
they lead to a smooth crossover between the dynami- 
cal states, with intermittent tumbling and tank-treading 
motions. Our simulations are in semi-quantitative agree- 
ment with a theoretical description based on a general- 
ized Kcllcr-Skalak approach. 

Second, we have investigated the behavior of vesicles 
near walls. Close to a wall, tumbling is strongly sup- 
pressed. Furthermore, the vesicle is repelled from the 
wall by the hydrodynamical lift force. We have found 
by boundary-integral calculations that the hydrodynamic 
lift force decays with increasing wall distance y CTa like 
1/ (2/cm lnj/cm)- However, for small wall distances - in 
particular in the regime of the MPC simulations - an 
effective dependence is observed. With increasing 
viscosity contrast, the lift force becomes weaker, as the 
vesicle becomes less deformable. The lift force also de- 
creases with increasing reduced area A* , and vanishes in 
the circular limit. We find that our numerical data are 
well described by a y/1 — A* dependence. 

Our results show that there is a different behavior of 
the lift force at intermediate and large distances from the 
wall, and that the lift force decreases significantly with 
increasing viscosity contrast. This may shed some light 
on the behavior in three dimensions, where experiments 
show a dependence of the lift force on the wall-membrane 
distance h, which decays as hr 1 for distances smaller 
than the vesicle radius [301 ] , whereas a y~£ decay has been 
found theoretically in a small range 1.1 < y C m/Rp ^ 1-25 



of wall distances (28|. Thus, we hope that our results 
will stimulate new experiments and simulations in 3D 
over a wider range of wall distances, reduced volumes, 
and viscosity contrasts. 
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APPENDIX A: DERIVATION OF 2D 
GENERALIZED KELLER-SKALAK THEORY 

1. Keller-Skalak Theory in Two Dimensions 

Keller and Skalak derived analytical expressions for 
the inclination angle 6 and the average angular velocity u> 
for 3D vesicles of fixed ellipsoidal shape, with (x/ai) 2 + 
(y/0,2) 2 + (z/a^) 2 = 1, based on the Jeffery theory [52l |. 

Although the KS theory is formulated for vesicles in 
three dimensions, it is straightforward to transfer it to 
two-dimensional systems by simply taking the limit — > 
oo. The resulting cylindrical three-dimensional geometry 
is equivalent to a 2D vesicle with the shape of an ellipse, 
(x/ai) 2 + [yja-i) 2 = 1 with ai > 02- Let S' be the frame 
which has its origin at the center of the ellipse, and the x' 
direction points into the direction of the long axis. Then 
the local velocity v' of an element of the tank-treading 
membrane is assumed to be 

=w(-— 4.— 2a)- (Al) 

in the frame S'. We define the auxiliary variables 

_ 1 - ap 2 _ 1 - ap 2 _ 1 + Qp 2 

1 + «d 8«d 2 

where qp = (0,1+0,2) /(ai— 0,2)- The balance of torques on 
the membrane and the assumption that the work done on 
the vesicles by the shear flow is dissipated in the interior 
of the vesicle leads to the non-linear differential equation 

f = |[-l + i%D,A)co S (20)] (A2) 

^ D ,A) = / {/ 1+ 1 + / ^_ 1) }. (A3) 

Furthermore, the average angular velocity u> is found to 
be 

u cos(2fl) 

7 2/ 1 {l + / 2 (A-l)} 1 > 
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FIG. 12: (Color online) Free energy F(a) of the two- 
dimensional vesicle relative to the prolate shape for A* = 0.7, 
k/1 = 50fc B T, fc sp = W 4 k B T/a 2 , and k A = 80k B T/a 4 . The 
solid (red) and dashed (blue) lines represent the simulation 
data and fitted curve, respectively. 



2. Shape Equation in Two Dimensions 

The vesicle shape is expanded in Fourier modes with 
polar angle <j> as r(</>) = RAe r {l + J2 u m exp(im0/v / 27r)}. 
Based on the Stokes approximation and perturbation the- 
ory, the dynamics of a quasi-circular vesicle is described 
by Hi 



du m 
dt 



?7outi?A 



(A5) 



where T m = |m|/2(A 
3/2 + a), and h = 



l)(m 2 - 1), E m = (m 2 - l)(m 2 - 
2ir/2(\ + 1). A Lagrange multi- 
plier a keeps the perimeter L p constant. Following the 
procedure for 3D [l5j |. we decompose u± 2 into ampli- 
tude and phase, u±2 = r cxp(^f2i9) , and replace the 
force 2nE m r by dF/dr. Then, Eq. is obtained with 
a = 3r/^+0(r 2 ). 

The free energy F(a) for the same simulation parame- 
ters calculated with a version of the generalized-ensemble 
Monte Carlo method Q (see Fig. [T2"]) : The vesicle con- 
formations are sampled under the uniform distribution 
of a, and then the canonical distribution is obtaine d by 
the reweighing. Instead of an interpolation [f| 0, Il5j |. 
we use fit functions here to obtain smooth functions. 



The force is fitted 
9 + 180a - 110a 2 - 
relation an = 2a/ 3 
A* = 0.7. 



as a function —(l/k-oF)dF/da = 
- exp(80a — 46). We obtained the 
+ 0.14a 4 by fitting for the ellipse of 



APPENDIX B: 



MEMBRANE COLLISIONS WITH 
THE SOLVENT 



For the modeling of an impermeable membrane, inte- 
rior and exterior solvent particles have to stay on the 
appropriate side of the membrane. Depending on the 
position of the MPC particle with respect to membrane 
bonds, either one or two membrane monomers partici- 
pate in a membrane collision. The new velocities of the 



m outside 




FIG. 13: Spatial regions for the membrane collision between 
MPC particles and membrane monomers. MPC particles lo- 
cated in the dark-gray regions perform a two-body collision 
with the monomers i or i+. MPC particles in the light-gray 
regions collide with both monomers. For further explanations, 
sec text. 



"coil particles (the MPC particle and the n co \ 
branc monomers) after the collision arc then 



1 mem- 



2 (v cm + u x r ljC ) - ' 



(Bl) 



where v cm is the center-of-mass velocity of the n co u-body 
system, u) its angular velocity, and rj jC are the particle 
positions relative to the center of mass of the rt co n-body 
system. This is a bounce-back collision for the relative ve- 
locities which conserves both the total translational and 
angular momenta. The spatial regions for the selection 
of colliding MPC particles and membrane monomers are 
illustrated in Fig. [13] 
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